Mon. Not. R. Astron. Soc. 000, [TV?? (2010) Printed 14 January 2013 (MN MgX style file v2.2) 



Accretion states of the Galactic microquasar 
GRS 1758-258 



o 

(N 



in 



6 



> 
on 
o 
o 
cn 

cn 
o 



is 



Roberto Soria 1 ' 2 *, Jess W. Broderick 3 , JingFang Hao 4 , Diana C. Hannikainen 5 , 
Missagh Mehdipour 1 , Katja Pottschmidt 6 ' 7 and Shuang-Nan Zhang 8 

1 Milliard Space Science Laboratory, University College London, Holmbury St Mary, Surrey RH5 6NT, UK 

2 Gurtin Institute of Radio Astronomy, Curtin University, GPO Box U1987, Perth, WA 6845, Australia 
3 School of Physics & Astronomy, University of Southampton, Southampton, Hampshire SOU 1BJ, UK 
^Department of Physics and Tsinghua Center for Astrophysics, Tsinghua University, Beijing 100084, China 
5 Aalto University Metsahovi Radio Observatory, Metsdhovintie 114, FI-02540 Kylmdld, Finland 

6 CRESST and NASA Goddard Space Flight Center, Astrophysics Science Division, Code 661, Greenbelt, MD 20771, USA 

7 Center for Space Science and Technology, University of Maryland, 1000 Hilltop Circle, Baltimore, MD 21250, USA 

8 Key Laboratory of Particle Astrophysics, Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, China 



Accepted 2011 March 15 — Received 2011 March 05 — in original form 2011 February 1 



ABSTRACT 

We present the results of a radio and X-ray study of the Galactic microquasar GRS 
1758—258, using unpublished archival data and new observations. We focus in partic- 
ular on the 2000-2002 state transitions, and on its more quiet behaviour in 2008-2009. 
Our spectral and timing analysis of the XMM-Newton data shows that the source was 
in the canonical intermediate, soft and hard states in 2000 September 19, 2001 March 
22 and 2002 September 28, respectively. We estimate the disk size, luminosity and 
temperature, which are consistent with a black hole mass ~ 10M Q . There is much 
overlap between the range of total X-ray luminosities (on average ~ 0.02LEdd) in the 
hard and soft states, and probably between the corresponding mass accretion rates; 
in fact, the hard state is often more luminous. The extended radio lobes seen in 1992 
and 1997 are still present in 2008-2009. The 5-GHz radio core flux density has shown 
variability between ~ 0.1-0.5 mJy over the last two decades. This firmly places GRS 
1758—258 in the radio-quiet sequence of Galactic black holes, in the radio/X-ray plane. 
We note that this dichotomy is similar to the dichotomy between the radio/X-ray se- 
quences of Seyfert and radio galaxies. We propose that the different radio efficiency 
of the two sequences is due to relativistic electron/positron jets in radio-loud black 
holes, and sub-relativistic, thermally dominated outflows in radio-quiet sources. 
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1 INTRODUCTION 

Understanding the connection between black hole (BH) ac- 
cretion and jet properties is a fundamental astrophysical 
problem. It is the key for the modelling of BH growth and 
feedback effects. Much progress has been achieved over the 
last two decades, thanks to combined X-ray/radio stud- 
ies of transient stellar-mass BHs, over timescales of a few 
months/years. Taken together, X-ray spectral and timing 
properties, and radio jet properties are the e mpirical crite- 
ria that identify physical "accretion states" (|Fender et al.l 
12004 iMcClintock fc Remillardlbood h 

When the mass accretion rate is below a few percent 
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of the classical Eddington rate, the accretion flow is hot- 
ter, geometrically thick, can more easily support a poloidal 
magnetic field and accelerate/collimate outflows. Radiative 
efficiency decreases and mechanical efficiency (jet power) in- 
cre ases. Advection-dominated flows may s et in (ADAF mod- 
els: iNaravan &Ty3|1994 lEsin et aUll997t) . due to inefficient 
cooling at low densities. Observationally, the hard state is 
characterized by flat-spectrum radio emission from the core 
(direct signature of a continuous, self-absorbed jet), by an 
X-ray spectrum dominated by a hard power-law (photon in- 
dex r ~ 1.7), and by high short-term variability in the power 
density spectrum. 

At accretion rates above a few percent of the Edding- 
ton rate, most stellar-mass BHs are seen to switch to the 
thermal-dominant (high/soft) state, characterised by an op- 
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tically thick, geometrically thin disk. This state is radiatively 
efficient and the jet is suppressed: these BHs are radio quiet. 
The X-ray emission comes mostly from the inner part of the 
disk, near the innermost stable circular orbit. The inner ra- 
dius and peak temperature of the disk depend directly on 
the accretion rate and BH mass. 

However, the physical mechanisms that cause BHs to 
switch between the hard and soft state are still poorly un- 
derstood, as are the intermediate states between them. The 
sample of transient BHs (mostly Galactic BHs with a low- 
mass donor star), used to define the "canonical" state be- 
haviour may be biased. The luminosity threshold for the col- 
lapse of the thick flow and the suppression of the jet changes 
not only between different sources, but also for individual 
sources between different outbursts, and between the ris- 
ing and declining phase of the state transition cycle, due to 
hysteresis — once it is formed, it takes time for the accretion 
disk to dissipate after the accretion rate has dropped. Fur- 
thermore, some BH s (particula rly among the ultraluminous 
X-ray source class: ISoriall201ll ) do not appear to switch to 
the thermal-dominant state even at high luminosities, and 
have an X-ray spectrum persistently dominated by a power- 
law-like component. The origin of the power law component 
itself is still unclear, with possible contributions from ther- 
mal Comptonization in a hot corona, and sy nchrotron self- 
Com pton emission from the base of the jet l|Markoff et alj 
l2005h . 

In this paper, we investigate accretion states by study- 
ing the X-ray and radio properties of the Galactic mi- 
croquasar GRS 1758-258 (|Mandroul Il99d : ISunvaev et al.l 
ll99lD . which stands out among Galactic BHs for its persis- 
tent activity, mostly in a moderately luminous hard state. 
We focus on its behaviour during 2000-2002 (Figure 1), 
when it underwent a series of spectral state transitions, 
and in 2008-2009 (Figure 2), when it was in a more steady 
hard state. For the X-ray study, we use unpublished XMM- 
Newton data supplemented by Swift data. For the radio 
study, we use our recent observation from the Australia Tele- 
scope Compact Array (ATCA), supplemented by archival 
Very Large Array (VLA) data. 



2 MAIN PROPERTIES OF GRS 1758 258 

The BH candidate GRS 1758-258 is located in the direction 
of the Galactic centre, « 5° away from it on the sky. Here 
we shall assume a distance of 8 kpc. Its BH mass and orbital 
period are still unknown, because optical st udies of the stel- 
lar donor are hampered by high extinction dRothstein et alj 
l2002i : ISmith et alj|2002bl : iMunoz-Arionilla et alJl201Ch . The 
optical counterpart has a dereddened luminosity consis- 
tent wit h an early A-type main sequ ence star, but unusual 
colours (IMufioz-Arionilla et al.ll2010t) ; part of the optical/IR 
emission may also come from the accretion disk or the jet. 

GRS J 1758— 258 has a point-like core with radio and 
X-ray emission, and an extended lobe structure, with a pro- 
jected distance of w 2.5 arcmin between the lobes (sw 6 pc 
at d = 8 kpc). The radio core is variable or transient. For 
example, we already knew t hat it was undetected (< 0.1 
mjy) at 6 cm in early 1992 dRodrfguez et al.lll992l ). when 
the BH was in an X-ray soft state. It was detected at vari- 
able flux levels « 0.1-0.5 mjy in 1992 September, when the 
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Figure 1. Top panel: RXTE/PCA lightcurve of GRS 1758-258 
between 1999 and 2010, from the Galactic Bulge Scans, 
and timeline of our XMM-Newton and ATCA observa- 
tions. Bottom panel: zoomed-in view around the 2000-2002 
XMM-Newton observations, which coincided with a tran- 
sition to the soft state and back to the hard state. We 
have also marked the epochs when a compact jet was or 
was not detected in the VLA data. Galactic Bulge Scan 
lightcurves are provided by Craig Markwardt and are available at 
http: / /lheawww. gsfc.nasa.gov/users/craigm/galscan/main. html 



BH switched back to the h a rd state, in 1993 a nd in 1997 
dMirabel fc Rodriguez! Il994l : iMartf et al.l 120021 : lHardcastlel 
l2005h . The extended radio lobes, with an optically-thin syn- 
chrotron spectrum, appear to be persistent s tructures rather 
than fast-moving ejections (|Hardca stlc 2005). This is slightly 
unusual for Galactic microquasars ~(|Hemzll2002h . 

In contrast with a sparse radio coverage, there has 
been an almost continuous hard X-ray monitoring over 
the last two decades with GRAN AT CGRO, RXTE, I N- 
TEGRAL and Swift dLin et al J l200d: ISmith et all l2002al lbl: 



- _ j jft _ 

iPottschmidt et al.1 BoOrj ; lLohfink et al.l l20rih . There was 
also a series of four C handra observations in 2000-2001 
ijHeindl fc Smith! l2002ah. and three XMM-Newton obser- 
vatio ns in 2000-2002 jGoldwurm et all l200ll : iMiller et all 
2002). The X-ray data show that the source is always ac- 
tive (unusual for Galactic BHs with low-mass donor stars), 
with a variability by a factor of a few. It has spent most of 
the last two decades in a hard state, with occasional tran- 
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Figure 2. Swift/BAT lightcurve during 2008-2009, 
and timeline of the Swift/XRT spectral observations 
and of our ATCA observations. Swift/BAT lightcurves 
are provided by Hans Krimm and are available at 
http:/ /heasarc. nasa.gov/docs/swift/results/transients/ The 
scarcity or poor quality of X-ray coverage during most of Decem- 
ber and January is due to an unfavourable solar orientation. 



sitions to an intermediate or soft state. The typical 0.3-10 
keV luminosity is w 1-2 x 10 37 erg s _1 , corresponding to ac- 
cretion rates ~ 0.01-0.03 Eddington for typical stellar-mass 
BHs with masses ~ 10M© and radiative efficiency ~ 0.1. 



3 RADIO OBSERVATIONS 

3.1 VLA data from 2001-2003 

We searched the NRAO VLA Archive Survey database^ for 
observations of GRS 1758-258 during the 2000-2003 se- 
quence of state transitions. A total of 38 unpublished, pro- 
cessed radio maps at 4.86 GHz were found covering the pe- 
riod 2001 February 1 - 2003 November 21 (Table 1). The me- 
dian root-mean-square (rms) noise level near GRS 1758—258 
is about 75 /iJy beam -1 . We usually found that the best 
way to determine the core flux density was either to fit a 
point-source elliptical Gaussian, or to use the flux density 
of the peak pixel in the map. However, for the maps from 
2003 February 6 and 26, where it appears that the core is 
extended, we instead summed the pixels in a box or poly- 
gon surrounding the component. The angular resolution was 
generally too fine to obtain lobe flux measurements (the ex- 
tended emission was resolved out), and the few maps with 
lower angular resolution did not have enough sensitivity; 
therefore, we only used the VLA maps for studying the core 
emission. 

3.2 ATCA observations from 2008-2009 

We observed GRS 1758-258 at 6 cm with the ATCA in four 
separate sessions over the period 2008 December 29 - 2009 
January 02. The array configuration was 6C, with minimum 
and maximum baselines of 153 and 6000 m, respectively. We 

1 http: / /www. aoc.nrao.edu/~vlbacald/ 
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Table 1. Log of the VLA radio core observations at 4.9 GHz, 
between 2001-2003. The data for this table were obtained from 
images available in the NRAO VLA Archive Survey database. 
Upper limits are 3a levels. 



observed simultaneously in two bands centred at 4800 and 
5423 MHz, each with an effective bandwidth of 104 MHz. 
The total integration time on-source was 35.03 h, spread 
approximately evenly across the four sessions. The primary 
calibrator was B1934— 638, and B1752— 225 the secondary 
calibrator. 

Dat a reduction and imaging was carried out with 
miriad (|Sault et al.lll995h . We imaged all of the data to - 
gether using a robust weighting parameter |Briggsl Il995l ) 
of 0.5; the effective frequency of the combined data is 5.09 
GHz. In addition, the multi-freq uency deconvolution task 
mf 'clean jSault fc Wieringalll994h was used because of the 
frequency gap between the two bands. The cleaned, primary- 
beam-corrected image (Figure 3) has an angular resolution 
of 5'.'08 x l'.'69 (beam position angle 1?9). We also tapered 
the data to produce an image with the same resolution as 
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Table 2. Log of the XMM-Newton/EPIC+RGS observations 
(top) and Swift/XRT observations (bottom) used in this paper. 



the map made from the com bined 1992 and 1997 VLA data 
(Fig. 3 in lMarti et alll2002l ); the angular resolution in this 
case is 10'.' 1 x 6'.'0 (beam position angle 9? 7). A tapered map 
provides better sensitivity to the low-surface-brightness ex- 
tended emission from the lobes, which is resolved out in our 
higher-resolution image. In the vicinity of GRS 1758—258, 
the rms noise levels in the two maps are 16 jj,Jy beam - 
(higher resolution) and 26 /iJy beam -1 (tapered). More- 
over, we estimate that the calibration uncertainty relative 
to B1934— 638 is about 5 per cent in both cases. 



4 MAIN RESULTS OF OUR RADIO STUDY 

4.1 Radio core behaviour during 2001—2003 

The behaviour of the radio core during the 2001 soft state 
transition mirrored th at observed during the previous major 
soft transition in 1992 IjGilfanov et al.ll 19931 ; iGrebenev et al.l 
ll997l ; lKeck et al.ll200lj ). Between 1992 January and July, the 
core was undetected at 5 GH z, at a flux density < 0.07 mjy 
l|Mirabel fc Rodriguej[l994l ). It was d etected again on 1992 
Septe mber 10-26 at 0.53 ± 0.04 mjy (|Mirabel fc Rodriguez! 
1 1994! ) as the X-ray source returned to a harder state. It 
remai ned visible and varied betw een ~ 0.2-0.3 mjy during 
~ e3ll994l). 



1993 (iMirabel fc Rodriguez! |1994|). and between m 0.1-0.2 
mjy in 1997 (|Lin et alj|2000l ; iMarti et al.ll2002l ; lHardcastle! 
2005), while the X-ray source persisted in a hard state. 

Systematic VLA monitoring started again in 2001 
February (Table 1), when the core was not detected, down 
to a limit of 0.24 mjy (3 times the rms noise). At that epoch, 
the X-ray source was undergoing a series of rapid spec- 
tral state changes ([Smith et al.ll2001al ). The short transition 
timescale suggests that there was already a full disc, partly 
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Figure 3. Top panel: tapered radio map at 6 cm (5.09 GHz), 
from our ATCA observations on 2008 December 29-2009 Jan- 
uary 2. Contour levels are 65, 78, 91, 104, 107 fijy beam -1 , cor- 
responding to 2.5, 3, 3.5, 4, 4.5 times the rms noise. The beam 
size is lO'.'l X 6''0, and its position angle is 9?7. Bottom panel: 
higher-resolution zoomed-in view of the core region (source C) 
from the same ATCA data. Contour levels are 48, 80, 112, 144 
£tjy beam -1 , corresponding to 3, 5, 7, 9 times the rms noise. The 
beam size is 5'.' 09 X 1'.'69, and its p osition angle is 1?9 . Source D is 
proba bly a background source, see IMarti et akl j2002h ; [Hardcastlel 
l|2005h . 
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covered by a rapidly varying C omptonizing region (as dis- 
cussed in Sections 6.1, 7.2, and in lSmith et alj2~001af l. At the 
end of 2001 February, GRS 1758 -258 settled into a stead- 
ier soft state (|Smith et alj|200"lah . The radio core continued 
to be undetected during that month, down to < 0.1 mjy 
(Table 1). It was undetected throughout the rest of 2001, 
although the VLA observations were too short to place very 
stringent constraints. When the RXTE/PCA X-ray moni- 
toring resumed on 2002 February 5, the 2-10 keV count rate 
had already increased to a value close to its long-term aver- 
age (Figure 1), consistent with a return to an intermediate 
or hard X-ray state. However, the radio core was still unde- 
tected to < 0.24 mjy until at least 2002 March 30 (Table 1). 
The radio jet finally appeared in the 2002 July 12 snapshot, 
at a 5-GHz flux density m 0.27 mjy. It was then detected 
in every subsequent observation from 2002 and 2003 (peak 
pixel > 3.5cr in all cases), varying between sa 0.24-0.46 mjy 
(Table 1). 

Overall, the radio core behaviour is in agreement with 
the theoretical expectatio ns and with the beh a viour of other 
Galactic BH transients (| Corbel et al ] l200ll ; iFender et al.l 
l2004lBelloni et al.ll2010l ): the jet is suppressed when the X- 
ray source is in a disk-dominated thermal state. It appears 
that the jet had already vanished or at least strongly de- 
clined in the intermediate states immediately before and af- 
ter the soft state (2001 February and 2002 January-March). 
However, the sparse sampling and shallow VLA detection 
limits do not allow us to constrain the transition more 
strongly. 

4.2 Radio core and lobes in 2008-2009 

We found core emission centred at RA(J2000) = 
18 h 01 m 12 s 39i Dec (j2000) = -25°44'36'.'0 (Figure 1), with 
an error of 0'.'08 in RA and 0!.'2 in Dec. This is still consistent 
with the position determined from the highest -resolution 
VLA observations dMunoz-Arionilla et alj |2010| . and from 
our direct re-analysis of the VLA data). In the higher- 
resolution image, the integrated core flux density, deter- 
mined from an elliptical Gaussian fit with the miriad task 
imfit, is 0.25 ± 0.02 mjy. Similarly, the integrated core flux 
density in the tapered map is 0.25 ±0.03 mjy. Though there 
is a hint of source extension for the core component in the 
higher- resolution map (Figure 3, bottom panel), more sensi- 
tive data are required to confirm this. The location of the two 
extended lobes is also consistent with those mapped in the 
1997 observations (|Martf et al.ll2002l ; lHardcastle(l2005l ). The 
peak flux density in the northern lobe (also known as VLA 
source B) is 135^Jy beam" 1 at RA(J2000) = 18 h 01 m 13f:04, 
Dec(J2000) = -25°43'33'.'9. Its integrated flux density is 
(0.4 ± 0.1) mjy, which we determined by summing the pix- 
els in a box surrounding the lobe. The peak flux density 
in the southern lobe (VLA source A) is 103/iJy beam -1 
at RA(J2000) = 18 h 01 m lK37, Dec(J2000) = -25°45'50'.'8; 
we cannot get a reliable integrated flux density because the 
signal-to-noise level is too low. We do not have enough sen- 
sitivity in the tapered image to determine whether the lobes 
have changed in flux and morphology. The integrated flux 
density f or the northern lobe wa s listed as (0.52 ±0.05) mjy 
in 19 92 ilRodriguez et all I1992T ) and (0.65 ± 0.06) mjy in 
1997 |Marti' et alj|2002l ). We also note that the nearby VLA 
source D, located « 25" north-east of GRS 1758-258, well 



visib l e in 1992 and 1997 (iRodrfguez et al.lll992l :l Martf et al.l 
2002; lHardcastlel l2005h and interpreted as a background 
source, has faded somewhat; its peak flux density in the 
higher-resolution image is ~ 0.09 mjy beam -1 . 



5 X-RAY OBSERVATIONS 

5.1 XMM-Newton observations from 2000-2002 

GRS 1758-258 was observed with XMM-Newton on 2000 
September 19, 2001 March 22 and 2002 September 28 (Ta- 
ble 2): that is, immediately before, during and immediately 
after the 2001 soft state. An analysis of the 2001 Reflec- 
tion Grating Spec trograph (RGS) data was presented by 
iMiller et al.l l|2002l ). Here we present the results of the Eu- 
ropean Photon Imaging Camera (EPIC) observations in the 
three epochs, and of the 2002 RGS observation; we also re- 
examine the RGS spectrum from 2001 in view of the new 
EPIC results. 

We downloaded the Observation Data Files from the 
public archive and used the Science Analysis System (SAS) 
version 10.0.0 (xmmsas_20100423) to process and filter the 
event files and extract spectra and lightcurves. We selected 
single and double events (pattern ^ 4 for the pn and pat- 
tern ^ 12 for the MOS), with the standard flagging criteria 
#XMMEA_EP and #XMMEA_EM for the pn and MOS 
respectively. After building response and ancillary response 
files with the S AS tasks rmfgen and arfgen, we used XSPEC 
IjArnaudll 19961 ) version 12 for spectral fitting, and standard 
FTOOLS tasks for timing analysis. 

The observing modes for the three EPIC instruments 
reflect attempts to reduce a potential pile up problem. In 
2000, the pn camera was in small window mode, MOS1 in 
timing mode, and MOS2 in the non-standard prime partial 
refresh frame store (RFS) mode. In 2001, the pn was in large 
window mode, and both MOSs in partial window mode. In 
2002, the pn was in small window mode, MOS1 in timing 
mode, and MOS2 in partial window mode. The EPIC data 
are clearly affected by pile up: typical pn count rates in a 45" 
circle were w 100 ct s -1 and « 40 ct s -1 for the MOS in the 
first two observations, and about half that value in the third 
observations. We expect a deterioration in the response for 
count rates > 50 ct s -1 for the pn in small window mode, 
and > 6 ct s _1 in large window mode; pile up affects the 
MOS for rates > 5 ct s -1 in small window mode. The MOS1 
observations in timing mode from 2000 and 2002 do not 
suffer from pile-up. However, using timing data for spectral 
analysis introduces its own set of calibration uncertainties, 
that are probably as significant as the pile-up error. 

We used the SAS task epatplot to estimate the level of 
pile up for different source extraction regions, by comparing 
the observed and expected distribution of single and multi- 
ple events. After several tests, we chose to extract the source 
spectra from annuli between an inner radius of 10" and an 
outer radius of 45". We repeated the extraction and spec- 
tral analysis for source annuli with inner radii of 5", 15" 
and 20" and found that the difference in the fit parameters 
were within the error bars of the values for a 10" inner ra- 
dius. Thus, we cannot completely eliminate the errors due 
to pile- up, but we are confident that they are a small ef- 
fect compared with the spectral state transitions discussed 
in this paper. 
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Figure 4. Top panel: XMM-Newton/EPIC spectra from the 2000 
observation, fitted with an absorbed disk-blackbody plus power- 
law model; in the online colour figures, pn data are in red, MOS2 
in blue. The cyan and orange curves represent the power-law and 
disk-blackbody components in the best fitting pn model; see Ta- 
bles 3,4 for the fit parameters. Middle panel: XMM-Newton/EPIC 
spectra from the 200f observation, fitted with an absorbed disk- 
blackbody plus power-law model; pn data are in red, MOS2 in 
blue in the online version. The cyan and orange curves represent 
the power-law and disk-blackbody components in the best fitting 
pn model (Table 5). Bottom panel: XMM-Newton/EPIC spec- 
tra from the 2002 observation, fitted with an absorbed broken 
power-law model (Table 6); pn data are in red, MOS2 in blue in 
the online version. 
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Figure 5. Comparison between XMM- Newton/EPlC-pn spectra 
in 2000 (green in the online version; intermediate state), 200f 
(red; soft state) and 2002 (blue; hard state). 
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Figure 6. Top panel: XMM-Newton/RGS spectrum from 2001, 
fitted with an absorbed disk-blackbody plus power-law model. 
Bottom panel: XMM-Newton/RGS spectrum from 2002, fitted 
with an absorbed power-law model. See Table 7 for the best- 
fitting parameters. 
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Another calibration issue we had to deal with is the use 
of the RFS read-out mode for the 2000 MOS2 observation. 
This mode has been used only for few observations with high 
count rates (e.g., the Crab Nebula and the N132D supernova 
remnant) at very early stages of the XMM-Newton mission. 
The flux in this observing mode has never been properly 
calibrated. After consultations with the calibration scien- 
tists, and from a comparison with those early observations, 
we concluded that the RFS mode underestimates the true 
flux by about 30%, but the spectral shape is mostly reliable, 
consistent with spectra obtained in other modes. To double- 
check that, we fitted the 2000 MOS2 spectra twice: once 
independently, and once with fit parameters fixed equal to 
the pn parameters but with a free normalization constant 
(which turned out to be « 0.7, as expected). 

For the pn observation in large window mode, we ex- 
tracted the background from three 45" circles at sufficiently 
large distances from the source. For the 2000 and 2002 pn 
observations, in small window mode, we extracted the back- 
ground from two circular regions of 45" radius, located at 
the bottom of the observing window, ~ 3' from the source 
position. However, we checked that the background count 
rate is only modestly significant for the 2000 pn observation, 
when it is ~ 2% of the source count rate, and slightly affects 
the spectral fit only at energies > 10 keV. In 2001 and 2002, 
the background count rate is < 1% of the source rate, and 
its effect on the spectral shape is much less than the errors 
due to calibration uncertainties and pile-up corrections. For 
the MOS observations, we could not extract a background 
from the central chip because the whole observing window is 
affected by the source emission. From the count rate in the 
outer chips, we estimated that the MOS background rate 
in our fitting band is < 1% in all three epochs, and can be 
safely ignored. 

For all spectral fits, we rebinned the spectra to a min- 
imum of 30 counts per bin, so that we could use Gaussian 
statistics. For plotting purposes only, we rebinned the data 
to a minimum signal-to-noise ratio of 12 for the pn and 6 
for the MOS spectra. We fitted the pn spectra in the 0.8-12 
keV range, and the MOS spectra in the 0.8-10 keV range. 
We plotted the spectra down to 0.4 keV in Figure 4 for il- 
lustration purposes only. The small discrepancy between pn 
and MOS at 0.4-0.8 keV may be due to a number of reasons: 
a different level of pile-up, a higher contamination from ex- 
tended thermal emission or nearby soft sources in the pn 
(which has a larger point spread function), a more complex 
structure of the absorbing medium, or cross-calibration un- 
certainties. A soft excess below .9 keV is also clear ly visible 
in the RGS spectra from 2001 (|Miller et al.ll2002h . In any 
case, these are minor uncertainties, irrelevant to the main 
objectives of our spectral study. 

The XMM-Newton/KGS data were also processed and 
analysed from the Observation Data Files with the SAS 10.0 
software. Both RGS detectors were operated in the stan- 
dard Spectro+Q mode for the 2001 and 2002 observations. 
(Unfortunately, during the 2000 observation both RGS1 and 
RGS2 were inactive and thus no RGS data were obtained.) 
We filtered out time intervals with count rates > 0.3 ct s _1 
in CCD number 9, which is the most affected by background 
flares. For the 2001 observation, we removed « 3 ks, leaving 
a net exposure time of w 18.3 ks for each RGS detector. For 
the 2002 observation, the background level was lower and no 



time filtering was required; the net exposure time is m 8.5 
ks. After extracting the source and background spectra, we 
created response matrices with rgsproc. For the purpose of 
spectral fitting, we combined the first-order spectra and re- 
sponse matrices of RGS1 and RGS2 for each observation, 
using the SAS task rgscombvne. The channels in the RGS 
spectra were then rebinned by a factor of 2 to increase the 
number of counts per bin. We used the 0.8-2.1 keV energy 
range (10-26 A) for spectral fitting in xspec. 



5.2 Swift/XRT observations in 2008-2009 

We want to relate our ATCA radio observations of 2008- 
2009 to the X-ray spectral state at the time. To do that, we 
looked for archival X-ray observations as close as possible to 
that epoch. Swift/XRT observed the source four times be- 
tween 2008 April 01-11, and twice on 2009 March 12-13 (Ta- 
ble 2; see NASA's Heasarc data archive for a more detailed 
logbook). We used the on-line XRT data product generator 
(|Evans et all 120091 to extract lightcurves and spectra, in- 
cluding background and ancillary response files; we selected 
grade 0-12 events. We downloaded the suitable spectral re- 
sponse files for single and double events in photon-counting 
(PC) mode and window-timing (WT) mode from the latest 
Calibration Database (2009 December 1). 

In both the 2008 and 2009 observations, the XRT count 
rate is high enough (a few counts s~ ) to create serious pile- 
up problems in PC mode, but not in WT mode. To atten- 
uate the pile-up p roblem, the source extraction procedure 
l|Evans et alj 120090 first determines whether the count rate 
within a 30 pixel (70'. / 8) radius is > 0.6 counts s _1 : that 
is the level at which the observed single-pix el event rate is 
» 90 % of the input single-pixel event rate l|Vaughan et al.l 
2006). Then, a point-spread-function profile is obtained and 
compared with a default non-piled-up profile. From this, the 
extraction procedure determines the radius within which the 
two profiles differ significantly. Finally, the source extraction 
region is defined as the annulus between that radius and 
70'.' 8. 

The WT-mode spectra do not suffer from pile up (it 
would only be significant at > 100 counts s _1 ). However, 
the lack of spatial information does not allow background 
subtraction. That appears to be a problem at low energies 
(< 1 keV), for the particle induced background. In fact, at 
every epoch when both PC and WT spectra are available, 
the latter show a low-energy excess with respect to the PC, 
which is most likely due to soft proton contamination rather 
than true emission from the source. At energies > 1 keV, PC 
and WT spectra agree within the observational error. 

We examined the four individual spectra from the 2008 
and the two spectra from 2009, and found that the multiple 
spectra from each epoch were consistent with each other. 
The Swift/TSAT lightcurve (Figure 2) confirms that there 
were no unsual flares or state transitions during the two 
epochs. We coadded the spectra from each of the two epochs 
to increase the signal to noise. We rebinned and fitted the 
coadded spectra with XSPEC Version 12. 
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Figure 7. Power density spectra in 2000, 2001 and 2002. In each 
epoch, the red datapoints (shown in colour in the online ver- 
sion) are the power density spectra extracted from the EPIC pn 
lightcurves; the green datapoints are extracted from the MOS1 
data. The three spectra are typical of the intermediate state, 
high/soft state, and low/hard state respectively. 
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Table 3. Best-fitting spectral parameters for the 2000 September 
19 XMM-Newton/EPIC o bservation. Units: jVg .Qal (line-of-sight 
Galactic extinction, from Kalberla ct al. 2005) and Ng - lnt (in- 
trinsic extinction) are in units of 10 22 cm -2 ; fcT^bb in keV; N pa 
in 10 _1 photons keV -1 cm -2 s — 1 at 1 keV; the observed flux 
is in 10~ 10 erg cm -2 s — 1 ; the emitted luminosities in 10 37 erg 
s —1 . Fg 3_ 12 is the disk-blackbody fraction of emitted luminos- 
ity in the 0.3—12 keV band. Errors indicate the 90% confidence 
interval for each parameter of interest. As explained in the text 
(Section 5.1), the MOS2 flux values are 30% lower than the 
more reliable pn values. 



6 MAIN RESULTS OF OUR X-RAY STUDY 

6.1 X-ray spectral state in 2000-2002 

The main finding of our spectral analysis is that the source 
was in three different states during the three epochs: inter- 
mediate state (similar contributions from disk and power- 
law components) in 2000; soft state (disk dominated) in 
2001; and hard state (power-law dominated) in 2002. Here 
we will give a more detailed description of the spectral prop- 
erties in the three states. 

For the 2000 dataset, we used EPIC pn and MOS2 for 
the spectral analysis; we used MOS1 for the timing analy- 
sis (Section 6.2). We knew that the EPIC MOS2 calibration 
would underestimate the flux by « 30% (Section 5.1). We 
fitted the pn and MOS2 data in two different ways: first, with 
fully independent parameters, and then by locking the fit pa- 
rameters together but allowing for a free normalization fac- 
tor in the MOS2. We find that single component models (ei- 
ther power-law or disk-blackbody) are ruled out, but an ab- 
sorbed disk-blackbody plus power law model provides an ac- 
ceptable fit. Overall, the pn and MOS parameters are consis- 
tent with each other (Tables 3,4) , if we consider all the uncer- 
tainties outlined in Section 5.1. The power-law photon index 
r ~ 2.2; the disk colour temperature kT ln « 0.48 keV; the 
disk normalization iVdbb = (nn/dio) 2 cos 8 « 700. Follow- 
ing a conventional procedure, for fitting purposes, we split 
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Parameter pn Value 



MOS2 Value 



Model: const*phabs*phabs* (diskbb+pow) 
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MOS2 Value 



Model: phabs*phabs*pow 
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kT dhh 
A^dbb 

r 



1.00 



0.69 
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Table 4. Best-fitting spectral parameters for the 2000 September 
19 XMM-Newton/'EPIC observation, using a simple scaling con- 
stant between pn and MOS2 to account for the known calibration 
discrepancy (Section 5.1). Units are as in Table 3. Errors indicate 
the 90% confidence interval for each parameter of interest. 
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Table 5. Best-fitting spectral parameters for the 2001 March 22 
XMM-Newton/EPIC observation. Units are as in Table 3. Errors 
indicate the 90% confidence interval for each parameter of inter- 
est. 
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Table 6. Best-fitting spectral parameters for the 2002 Septem- 
ber 28 XMM-Newton/EPIC observation. Units are as in Table 3. 
Errors indicate the 90% confidence interval for each parameter of 
interest. 
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the total col umn density into a fix ed Galactic "line-of-sight" 
component (|Kalberla et a l. 2005)0 and a free additional in- 
trinsic component. The total column density ps 1.6 x 10 
cm -2 is similar to the values reported in previous X-ray 
studies. The unabsorbed luminosity in the 0.3-12 keV band 
is ~ 1.9 x 10 37 erg s _1 , and the direct disk emission con- 
tributes for « 30% in this band. 

In 2001, all three EPIC instruments were set in spec- 
tral modes. We fitted them independently and found that 
they give a consistent result (Table 5) . The X-ray spectrum 
is again well fitted with an absorbed disk-blackbody plus 
power-law model. This time, the disk-blackbody component 
dominates, and the spectrum is consistent with a canonical 
high/soft state. We are unable to remove the discrepancies 
between the three instruments in some of the fitting param- 
eters. However, we need to take into account that the pn 
(large window mode) was more heavily piled up than the 
MOS — hence, a larger fraction of its soft flux is incorrectly 
shifted to higher energies (cf. ^0.3-12 in Table 4). Further- 
more, MOS2 had several columns, straddling the source ex- 
traction region, that were nominally unflagged but returned 
an anomalously low count rate compared with other regions 
symmetrically placed with respect to the source position. 
Based on these considerations, we believe that MOS1 pro- 
vides the most accurate spectral fit parameters in 2001. In 
summary, the power-law index V ~ 2.9, the disk colour 
temperature kT[ n ~ 0.45 keV, and the disk normalization 
~ 1700. The direct thermal emission from the disk con- 
tributes m 80% of the 0.3-12 keV luminosity. Luckily, there 
is a VLA observation taken almost simultaneously with the 
XMM-Newton observation, on 2001 March 22: the radio core 
was undetected down to ~ 0.1 mjy (Table 1), which con- 
firms that the source was in the thermal-dominant state. 

In 2002, we used EPIC pn and MOS2 for our spectral 
analysis. The pile-up problem was somewhat less severe, and 
we found no other instrumental problems. We fitted the data 
simultaneously, keeping the fit parameters constant apart 
from a free normalization factor (Table 6). We found that a 
simple power-law with V ~ 1.5 is an acceptable fit. Adding a 
disk-blackbody component does not provide a significant im- 
provement to the fit, and the disk temperature remains very 
low and essentially unconstrained (also because of the high 
column density). We obtain a moderate, statistically signif- 
icant improvement (F statistic value = 29.4) by replacing 
the single power law with a broken power law (Table 6), 
with a slope r « 1.5 at E < 7 keV, and T w 1.9 at higher 
energies. Equally good fits are obtained using Comptoniza- 
tion models that account for the small degree of curvature 
in the spectrum. The closest VLA observation to this epoch 
was taken on 2002 September 18, and shows radio flux from 
the core, at » 0.24 mjy (Table 1); this confirms that the jet 
had started again. 

We tried to estimate a "total" X-ray luminosity by 
extrapolating the EPIC spectra to 0.3-100 keV in the 
three epochs. We have strong evid e nce from other satellite 
observations |Gilfanov et al.l 1 1995 ; iKuznetsov et aL 1 199 



iLin et all [20001 ; iKeck et all l200ll ; 



Pottschmidt et al 



2006) 



2 Although in this case we are only looking half way through the 
Galaxy, so the true line-of-sight extinction may be about half of 
that value. 



that the power-law spectrum extends to > 100 keV, at least 
in th e hard state; on on e occasion, it extended to » 500 
keV (|Lohfink et alj|201lf ). The 0.3-100 keV luminosity was 
w 2.3 x 10 37 erg s" 1 in the intermediate state, w 1.2 x 10 37 
erg s _1 in the soft state, and ~ 1.6 x 10 37 erg s _1 in the 
hard state. In the hard state, part of the accretion power is 
transferred as kinetic energy to the electrons in the corona 
and/or jet. So, hard states should be less radiatively effi- 
cient than the thermal-dominant state. Given the inferred 
luminosities, it is very likely that the mass accretion rate 
M was also slightly higher in the hard state than in the soft 
state (unless the former was moderately beamed, see Section 
7.2). This is unusual for Galactic BHs, as pr eviously noted 
llSmith et al. I l2001al ; IPottschmidt et"al]|2006h . Also unusual 
is the fact that over 20 years, the total X-ray luminosity has 
varied only by a factor of 3. 

In all three epochs, the pn and MOS1 timing-mode spec- 
tra have a small excess at energies < 0.8 keV, but that is 
not seen in the other MOS spectra, which have a higher spa- 
tial resolution and lower background contamination. This 
residual emission is inconsistent with the highly absorbed 
blackbody plus power-law components used to model the 
spectrum at energies > 0.8 keV. It can be formally mod- 
elled with an additional optically-thin thermal plasma com- 
ponent at kT < 60 eV behind a line-of-sight column den- 
sity. We suggest that these residuals are due to imperfectly- 
subtracted diffuse background emission, and arc unrelated 
to the point-like emission from GRS 1758—258, but we do 
not have enough data to prove it, also taking into account 
the additional uncertainties due to pile-up. 

The RGS spectra (Figure 6) are consistent with the in- 
terpretation resulting from the EPIC spectra. In 2001, the 
RGS spectrum is dominated by a thermal component with 
a small power-law excess at high energies. The RGS energy 
band is not sufficient to constrain the power-law component, 
so we fixed its slope and normalization to the best-fitting 
values derived from the MOS1 data (Table 5). We obtain 
that the best-fitting inner-disk temperature, normalization 
and total column density are also the same as in the MOS1 
spectrum (Table 7). The sma l l diffe rence between our re- 
sults and those of iMiller et ail (|2002D is due to our different 
choice of the power-law component. In 2002, there is no need 
for a thermal component: the RGS spectrum is well fitted 
with an absorbed power law (Table 7). The photon index of 
the RGS spectrum is slightly steeper than in the simultane- 
ous EPIC spectra, but marginally consistent considering the 
large error. Hints of a small soft excess below 0.9 keV (both 
in 2001 and 2002) are consistent with what is seen in the 
EPIC-pn spectra. We do not detect any significant features 
emission lines either in 2001 or in 2002. 



6.2 X-ray time variability in 2000-2002 

We examined the variability of GRS 1758—258 in the three 
epochs. For each epoch, we used MOS1 (in timing mode for 
the 2000 and 2002 observations) and pn (which provides the 
highest count rate). We extracted lightcurves binned to 0.01 
s, from the same regions used for the spectral analysis. We 
computed the Fourier transforms (powspec task in FTOOLS), 
normalized the resulting power spectra to be the fractional 
mean squared variability per unit frequency interval, and 
subtracted the Poisson noise level. 
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Parameter 


2001 Value 


2002 Value 


N B ,Gel 


0.75 


0.75 


N H ,int 


1 00+ 004 
1 - uu -0.04 


1 04 +0 18 


kT dhh 


°- 45 -0.02 




N Ahh 






V 


[2.85] 


2 02+ ' 47 


N po 


[0.054] 


1 7 +o.7 

i -'-0.5 



7Cu 



1.24 (580.0/469) 1.02 (477.7/469) 



Table 7. Best-fitting spectral parameters for the 2001 and 2002 
XMM-Newton/RGS observations. For the 2001 spectrum, the 
slope and normalization of the power-law component were fixed 
at the values found from MOS2. Units are as in Table 3. Errors 
indicate the 90% confidence interval for each parameter of inter- 
est. 



In the 2000 dataset, we found significant short-term 
variability below 10 Hz, with a power-law red noise com- 
ponent becoming dominant below ~ 10 -3 Hz, a fiat top 
between ~ 10 _3 -1 Hz, and a steepening above ~ 1 Hz. We 
do not have enough signal to noise to determine whether 
th ere are qu a si-per iodic oscillations at ~ 1-5 Hz, as found 
bv lLin et af] l|2000l ) from the 1997 RXTE observations, but 
the MOS1 power spectrum hints at their presence (Figure 7, 
top panel). The rms fractional variability between 10 -4 and 
50 Hz is (10 ± 1)%. In 2001, the short term variability was 
much reduced: the rms fractional variability is consistent 
with zero. There may be a weak power-law noise component 
similar to the one seen in the 2000 lightcurves. The 2002 
observations show the highest degree of short-term variabil- 
ity, with rms = (28 ± 1)%. The strong band-limited noise 
has a break at « 0.4 Hz. In summary, the timing properties 
of GRS 1758—258 in the three epochs are consistent with 
those of the canonical intermediate state (probably soft in- 
termediate), high/soft state and low/hard state respectively 
llBelloni et all Il996l; iMendez et all Il998l ; iReig et all 120031 : 
iMotta et al.ll2009l : iBelloni et alJ^Oia vT in perfect agreement 
with the spectral identification. 



6.3 X-ray spectral state in 2008 2009 

The Swift/BAT hard X-ray monitoring during 2008-2009 
(Figure 2) suggests that the source was in the hard state 
during that time, before and after our ATCA observation. 
The discontinuity in the X-ray monitoring between Decem- 
ber and January is due to the unfavourable solar orientation, 
and does not represent significant evidence of state transi- 
tions. The X-ray hard state is consistent with our radio core 
detection between the two epochs (Section 4.2). 

The coadded 2008 April spectra have an on-source time 
of 21.1 ks in PC mode and 1.7 ks in WT mode. Both spectra 
are well fitted above 0.9 keV by a simple power-law model 
with photon index P w 1.5 (Table 8, Figure 8) and total 
column density ~ 1.5 x 10 22 cm -2 . In fact, a simple power 
law is a good fit for all the (background-subtracted) PC- 
mode data, down to 0.5 keV. In WT mode, there is a small 



excess at energies < 0.9 keV, similar to the excess found in 
some XMM-Newton spectra. We attribute this excess either 
to unsubtracted particle background or to extended thermal- 
plasma emission around the source. Overall, the spectral 
parameters are very similar to those found in 2002 (Table 6 
and section 6.1), but the luminosity is a factor of 2 lower. 

The 2009 March spectra were taken only in WT mode. 
Considering the systematic uncertainty of using a timing 
mode for spectral analysis, we found that both a simple 
power-law and a broken power-law model provide acceptable 
fits to the data above 0.9 keV (Table 9, Figure 8); the broken 
power-law is preferred, suggesting mild spectral curvature. 
For this reason, equally good fits can be obtained either with 
phenomenological combinations of a straight power-law and 
broad-band curved components in the XSP EC arsenal, or 
with (more physical) Comptonization models (|Lohfink et al.l 
I2011J) . In any case, the photon index Y < 2 over the whole 
fitting band is again typical of the hard state. There is, 
again, a small excess at energies < 0.9 keV. That soft com- 
ponent is not consistent with a disk-blackbody spectrum 
behind the same column density as the power-law emission. 
As noted before (Section 6.1), it can formally be fitted with 
an optically-thin thermal plasma component at fcT < 60 eV 
and only line-of-sight column density. We do not think that 
this soft excess comes directly from GRS 1758—258. 

In both the 2008 and 2009 Swift datasets, the total X- 
ray luminosity was « 10 37 erg s -1 , a factor of 2 lower than 
in the 2002 hard state observation. Therefore, in this case, 
the hard state luminosity was slightly lower than the 2001 
soft state luminosity, unlike what is seen in other (more lu- 
minous) hard state observations of this BH. This confirms 
that the spectral state of this BH is not solely determined by 
the total luminosity or mass accretion rate. Based on these 
data alone, we cannot speculate about the physical struc- 
ture of the accretion flow in the 2008-2009 hard state: e.g., 
whether it was a truncated disk plus advective region, or a 
full disk covered by a hot corona or outflow. T his issue will 
be investigated elsewhere (|Lohfink et al.ll201ll ). 



7 DISCUSSION 

7.1 Disk size and BH mass 

The peak colour temperature of the disk is similar (fcTi n ~ 
0.45-0.48 keV) in the two epochs of our XMM-Newton ob- 
servations when a disk component is directly visible (Tables 
3,4,5). A colour tempera ture 0.48 keV was also found by 
IPottschmidt et alJ {2006) from INTEGRAL observations in 
2003 March-April, another epoch when the X-ray spectrum 
was dominated by the soft thermal component. RXTE/PCA 
observations from 2001 February-March (a few days before 
our 2001 observations) also s howed a dominant d isk compo- 
nent with fcTi n ~ 0.46 keV (|Smith et al.ll2001al ). A Chan- 
dra/ ACIS HETG spectrum from 2001 M arch 24 showed 
a peak colour temperature ~ 0.51 keV (jHeindl fc Smith] 
l2002bl ). 

The presence of a soft excess was suggested even when 
the source was in a harder, power-law dominated state. 
ROSAT/PSPC observations from 1993 March 31-April 1 are 
consistent with a thermal component wi th fcTbb ~ 0.5 keV 
|Mereghetti et al.lll994l : lKeck et al.ll200ll ): in the ASCA ob- 
servations of 1995 March 29-30, there are hints of a heavily 
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Figure 8. Top panel: coadded S«/i/f/XRT spectra from the four 
observations between 2008 April 01—11. The spectrum in PC 
mode is plotted in red in the online version; the WT mode in 
blue. The model is an absorbed power-law, fitted between 0.9-10 
keV for both spectra. See Table 8 for the best fitting parame- 
ters, and Section 6.3 for an explanation of the soft excess in the 
WT spectrum. Bottom panel: coadded Sm/i/XRT spectrum (WT 
mode) from the two observations of 2009 March 12-13. The model 
is a broken power law, fitted between 0.9-10 keV. See Table 9 for 
the best fitting parameters. The soft excess is formally well fit- 
ted by a low-temperature thermal plasma component with only 
line-of-sight absorption (plotted in pink in the online version); see 
Section 6.3. 



Com ptonized seed thermal c omponent with fclbb ~ 0.4-0.5 
keV (fMcreghe tti et aD 1 19971 ). However, both observations 
can also be fitted with simple power laws ijMereehetti et aD 
ll997l ; lKeck et al.ll200lh . 

The simplest explanation for this consistent value of 
the temperature is that it corresponds to the peak colour 
temperature of a full disk, extending all the way to the in- 
nermost stable circular orbit. If so, Ti n ~ L 1 / 4 . Since the 
X-ray luminosity and probably the total accretion rate have 
varied only by a factor < 3 over the last two decades, the 
peak temperature of a full disk (not significantly hardened 
by Comptonization) should vary only by a factor < 3 ' 25 , 
that is only < 30%. Most Galactic BH transients have typ- 
ical disk temperatures ~ 1 k eV in their high/soft state 
jRemillard fc McClintock1l2006h . with X-ray luminosities ~ 



Parameter 


PC Value WT Value 


Model: const*phabs*phabs*pow 


Const 

W H ,Gal 

r 

JVpo 


1.00 1.15«;° 4 
0.75 
0.681°;°? 

49+ 04 

—0.04 


/o.3-12 
£n. 3-12 
Lo. 3-100 


3 1+ - 1 3 6+ ' 1 
a - L —o.i °' u -0.1 

35+ ' 01 41+ - 01 
u * ao -o.oi u -^ ± -o.oi 

-, -,+0.1 -i ,+0.1 
1 ' 1 -0.1 1,o -0.1 


X u 


0.90 (614.1/683) 



Table 8. Best-fitting spectral parameters for the coadded 2008 
April SW/t/XRT observations. We allowed for a free normaliza- 
tion parameter between the PC and WT spectra. Units are as 
in Table 3. Errors indicate the 90% confidence interval for each 
parameter of interest. 

a few 10 38 erg s" 1 . GRS 1758-258 has an X-ray luminos- 
ity ~ 2 x 10 37 erg s _1 ; if its BH has a similar mass to the 
other Galactic BHs, we do indeed expect a disk tempera- 
ture f» 0.5 keV when all or most of the accretion power is 
directly radiated by the disk. That is the only state when 
we can constrain the BH mass from the disk temperature. 

We can use the disk normalization parameter iVdbb for 
a more quantitative estimate of its BH mass. The best-fit 
value iVdbb « 700 found in the 2000 XMM-Newton dataset 
certainly underestimates the thermal emission from the full 
disk, because at least half of the disk photons have been 
upscattered into the power-law. This would not change the 
observed disk temperature but would reduce its flux nor- 
malization. In 2001, the disk component dominates and our 
most reliable spectral fits suggest iVdbb ~ 1700. The 2003 
INTE GRAL study gives AT dbb « 2700 during another soft 
state l|Pottschmidt et afll2006l ). Let us take an average value 
of 2200 for the sake of our ballpark mass estimate: then, the 
apparent inner radius r- m = riio(A?dbb/ cos 6*) ' 5 ~ (40/ cosf?) 
km. The physical inner radius R ln w 1.19ri n ~ (45/ cos 9) 
km (|Kubota et al.|[T998l ). The inclinaton angle 9 is unknown, 
but the source is unlikely to be perfectly face-on (given 
the extended jet lobe structure), or edge-on (because no X- 
ray eclipses have been seen). Assuming a Schwarzschild BH 
(Rin ~ 6M) we obtain a characteristic mass ~ (5/ cos 9) Mq, 
typical of Galactic stellar-mass BHs, and therefore an aver- 
age X-ray luminosity ~ 2% Eddington. 

7.2 Disk and corona during state transitions 

Taken individually, all three main spectral states of GRS 
1758—25 8 are similar to the canonical states of Galac - 
tic BHs (|Remillard fc McClintock1l2006l ; iFender et al.ll2004r ). 
The characteristic long-term X-ray luminosity w 2 x 10 37 erg 
s _1 is also consistent with the characteristic threshold for 
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Parameter 


Value 


Model: phabs*phabs*pow 


N H dnt 

r 


0.75 

71+ - 06 

1 59+ 004 

69+ 04 
u - oa _ 0.04 


/0.3-12 
^0.3-12 
^0. 3-100 


„ 7 +0.1 
J - '-0.1 

45+ 01 

U -^° — 0.01 

1 2+ ' 1 


Xv 


1.18 (619.3/527) 


Model: phabs*phabs*bknpow 


N H dnt 

r 2 


0.75 

36+ 10 

1 12+ - 11 

3-8t°o:! 

1.88±S;S 
40+ 06 


/0.3-12 
^0.3-12 
Lo. 3-100 


o c+0.1 

J -°-o.i 


xl 


1.02 (536.8/525) 



that is often taken as the standard paradigm for hard/soft 
transitions in Galactic BHs. 



Table 9. Best-fitting spectral parameters for the coadded 2009 
March Swift/XRT observations. Units are as in Table 3. Errors 
indicate the 90% confidence interval for each parameter of inter- 
est. 



state transitions. However, the transition betw een states in 
GRS 1758—258 is unusual, as pre viously noted (|Smith et al.l 
l2001al : iPottschmidt et al.1 120061) . The radiatively-efficient 
soft state has a slightly lower total luminositj0, and hence 
lower mass accretion rate, than during most (not all) of the 
long-duration hard state. For example, our XMM-Newton 
spectral analysis shows that during the 2001 soft state, the 
source was slightly less luminous than in the harder states 
of 2000 and 2002, although slightly more luminous than 
in the 2008-2009 hard state. Moreover, the disk compo- 
nent has been seen to appear and disappear over very short 
timescales (Table 10), but always with a peak colour tem- 
perature « 0.45-0.50 keV. This is d ifficult to reconci l e with 
the disk truncation/refill scenario (|Esin et al.l Il997l , [l998) 



3 The dramatic drop in the 2-10 keV flux detected by RXTE 
(Figure 1) is largely compensated by the increase in flux from the 
disk, below 2 keV. 



As an alternative explanation, ISmith et al.l (|2001al . 

2002a) suggested that a full accretion disk is already present 
before the hard-to-soft transition, and is still present at least 
for some time after the return to the hard state. In this sce- 
nario, sometimes the disk is directly visible, but more often 
it is completely covered by a Comptonizing medium above 
it. The presence of full disks with inner radii close to the in- 
nermost stable circular orbit even in the low/har d state, at 
lumin osities ~ 0.01 Eddington, was suggested bv lReis et al.l 
(2010) from their study of several Galactic BH transients. 
Their findings suggest that transitions from the high/soft 
state to the low/hard state may be driven by changes in the 
corona (perhaps related to jet formation) rather than by the 
disappearance of the accretion disc. The behaviour of GRS 
1758—258 suggests the same scenario for transitions in the 
both directions, from hard to soft and vice versa, over short 
periods of time. 

One mechanism to produce a rapidly chan ging corona 
over a steady disk is the two-flow scenario l|Smith et al.l 
l2001al . | 2002a; Chakrabarti & Titarchuk 1995): a keplerian 
accretion flow feeds the disk (evolving over a viscous 
timescale) and an independent sub-keplerian flow feeds the 
hot Compton cloud (evolving over a dynamical timescale). 
Soft states correspond to the temporary shutoff of the sub- 
keplerian flow. One difficulty of this scenario is that the 
donor is a low-mass star, an unlikely source for a low- 
angular-momentum component of the accretion flow. 

We suggest another possibility: the hard emission com- 
ponent may come from a magnetically-powered coronal out- 
flow. Radio studies provide independent evidence that there 
is a jet during the hard states. Coronal magnetic fields may 
dissipate a large fraction of accretion power and transport 
it out as Poynting flux or as a magnetically dri ven wind 
|Merloni fc Fabianl [20021 : iKuncic fc Bicknelll r2004r i . In this 
scenario, the accretion flow may switch between a stan- 
dard disk and a coronal outflow in response to changes in 
the configuration of the poloidal magnetic field, and hence 
in the fraction of angular momentum transported by the 
azimuthal- vertical component of the magnetic stress rather 
than by the radial component. The underlying full disk 
may always be present even during the outflow dominated 
phases, but it will be masked by the outflow, and also colder 
and less luminous than in the disk dominated p hases, for 
the same mass inflow rate l|jollev fc Kuncicl l2008), Switches 
between a colder and a hotter disk phase may occur on 
the thermal timescale, which is shorter than the viscous 
timescale required to regrow a disk fro m large radii (i t h 



{H/R) 2 t viBC < i visc : iFrank et all 120021). Magnetohydrody- 



lagnctonyaroay- 

namic simulations sug gest (|De Villiers et al.l 120031 . 120051 : 
lHawlev fc Krolikl [2006) that magnetized outflows are not 
isotropic: they are confined by funnel walls. In this scenario, 
we speculate that for GRS 1758—258 we are looking down 
the funnel, so that the apparent luminosity of the hard state 
is enhanced by a factor of a few compared with the lumi- 
nosity in the soft thermal state (which has a more isotropic 
emission). This would make the true luminosity and accre- 
tion rate of the hard state lower than in the soft state, as 
observed in most other Galactic BHs. 
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Epoch 


Spectral State 


Disk Temperature 


Power-law Slope 


Radio Core? 


References 


2000 Jan-Aug 


Hard 


disc not seen 


S3 1.5-1.8 


? 


[1] 


2000 Sep 


Int 


S3 0.5 keV 


2.0-2.6 


? 


[1,2] 


2000 Oct-Dec 


Hard 


disc not seen 


S3 1.5-1.8 


? 


[1] 


2001 Jan 23 


Int 


S3 0.5 keV 


f» 2.2 


? 


[1.3] 


2001 Jan 29-Feb 21 


Hard 


disc not seen 


S3 1.8-2.0 


NO 


[1.3] 


2001 Feb 27 


Soft 


S3 0.5 keV 


S3 2.6 


? 


[3,4] 


2001 Mar 2-24 


Soft 


S3 0.5 keV 


S3 2.8-3.2 


NO 


[2,3,4,5] 


2001 Apr-Nov 


Soft 


? 


? 


NO 


[2,6,7] 


2002 Jan-Mar 


Int? 


? 


? 


NO 


[2,4,7] 


2002 Jul-2003 Feb 


Hard 


? 


? 


YES 


[2,6,7] 


2002 Sep 28 


Hard 


disc not seen 


S3 1.5-1.7 


YES? 


[2] 


2003 Mar-Apr 


Soft/Int 


S3 0.5 keV 


S3 2.3 


? 


[6,8] 



Table 10. Timeline of the 2000-early 2003 state tran sitions. "Disc not see n" means that a ther mal disk component i s not r equired to fit 



the X-r ay spec trum at those epochs. R eferen ces: [1] = [ Sm ith et a ll ll2002ah , [2 ] = this work, [3] =ISmith et al 
ll2001bl), [51 = iHeindl fc Smith! |2002b), [61 =|Pottschmidt et al.l <2006h. [71 = Heindl 



12] 

Smith (2002a), [81 



. fcOOlaT), [41 =ISmith et al.l 
Pottschmidt et al.l d2008f). 



7.3 Radio louder and radio quieter microquasars 

There is now strong evidence, from both the 1992 and 2001 
state transitions, that the jet in GRS 1758—258 switches off 
during soft states, similarly to most other stellar-mass BH 
transients. This is the basis of the canonical model for the 
disc-j et coupling in black hole X-ray binaries (jFender et al.l 
120041 ) and AGN (iKording et al.ll2006t ). When accreting BHs 
are in the jet-dominated state, there is a fundamental cor- 
relation between their X-ray luminosity and radio core lu- 
mino sity ([Merloni et al.ll2003l ; iGallo et a l. 2003; Falcke et al.l 
120041 ). As more and more Galactic BHs have been studied in 
radio and X-ray bands, it was found (|Ga llo 2007; Soleri et "all 
l20ld : ICalvelo et aT1l2010l ; ISoleri fe Fenderll201ll ) that there 
are two distinct groups or sequences of sources, with similar 
X-ray luminosity and hard spectral state but radio luminosi- 
ties differing by a factor ~ 10-30. 

GRS 1758—258 belongs to the radio-quieter group, de- 
spite its characteristic large-scale radio jet and lobes. In fact, 
this source is the radio quietest of all Galactic BHs in that X- 
ray luminosity range, S3 50-100 times fainter that the radio- 
loud sequence (Figure 9). Even if its X-ray luminosity were 
revised downwards by a factor of a few due to geometrical 
beaming (Section 7.2), GRS 1758-258 would still be in the 
radio-quieter group. 

The existence of radio-loud and radio-quiet stellar- 
mass BHs is reminiscent of the AGN classification. For 
a given range of black hole masses and accretion rates, 
som e kinds of active galaxies have stronger je ts than oth- 
ers (| Wilson fc Colbert)il995l ; ISikora et al.ll2007ft . In the case 
of AGN, this dichotomy corresponds to morphological dif- 
ferences in the host galaxy: elliptical galaxies are generally 
more radio-loud than disk galaxies. It is still not clear what 
corresponding environmental parameter is associated to the 
radio dichotomy in stellar-mass BHs: no obvious dependence 
was found on orbital periods, disc sizes, outburst history 
l|Soleri fc Fenderll201l| y 

Spin is often suggested as the intrinsic physical differ- 
ence between radio loud and radio quiet BHs: fast-spinnin g 
black holes can launch more powerful jets (|Blandfordllf990h . 
This is consistent with the scenario that black holes in gi- 
ant elliptical galaxies evolved to higher spins than those in 



spiral galaxies. However, we do not have enough reliable 
spin measurements for stellar-mass BHs to test whether 
the same correlation applies there. It is also still un- 
clear whether the jet power has a direct contribution from 
the BH spin jBlandford fe Znaiekl [19771: I McNamara et al 



l20ld:lFender et al.H2010l ). Furthermore, IBroderick fc Fendei 
|201lf ~ 



1 2011) suggest that BH spin has considerably less impact on 
AGN jet power than has been previously reported. 

Alternatively, the radio dichotomy in stellar-mass 
BHs may reflect di fferent modes of disc/jet coupling 
l|Soleri fc Fenderll201ll ). or different structures of the accre- 
tion flow in the hard state. For example, the efficiency of jet 
launching may be different between an ADAF-like configu- 
ration (truncated disk and hot, advective inner region) and 
a disc-plus-corona structure (corona on top of a full disk). 
Both accretion states may be phenomenologically classified 
as "low/hard state" with similar X-ray properties, but dif- 
ferent radio properties. 

We note that the observed presence of two distinct se- 
quences of Galactic BHs in the radio/X-ray luminosity plane 
(Figure 9) is similar to the dichotomy between radio galaxies 
(including both FRII and FRI objects) and Seyfert galax- 
ies. Both radio galaxies and Seyferts show a strong corre- 
lation between core radio and X-ray luminosity (Fig. 2 in 
iPanessa et al. 2007), or core radio and [OIII] luminosity 
(Fig. 4 in lBicknellll2002l ). which is another proxy for the 
radiative accretion luminosity. But the two sequences are 
shifted by almost three orders of magnitude, with Seyfert 
galaxies being fainter in the radio for a given X-ray luminos- 
ity. Even after both samples are normalized b y their masses, 
accor ding to the fundamental plane relations (| Merloni et al.l 
2003), they still form two distinct sequences, shifted by a 
factor ~ 50 in rad io luminosi ty (Fig. 7, right panel, in 
IPanessa et aT]|2007l ). Similarly, IBroderick fc Fenderl (|201ll ) 
plotted radio loudness (also determined using core radio lu- 
minosities, and with a mass correction applied) as a function 
of the Eddington ratio. They found that radio galaxies are 
on average more radio-loud than Seyfert galaxies and LIN- 
ERS by a factor of ~ 30 for Eddin gton ratios < 1 per cent. 

By analogy with the model of iBicknelil (|2002T ). we sug- 
gest that this dichotomy corresponds to the different radio 
efficiency of relativistic electron-positron jets (the radio loud 
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Figure 9. Radio loud and radio quiet sequences of Galactic BHs. 
The red box marks the characteristic range of X-ray and radio 
variability of GRS 1758—258 during the hard state (when the 
radio core is detecte d) . Datapoints from all other sources are from 
ICalvelo et all <201Ch . 



sequence) and sub-relativistic, thermally dominated bary- 
onic jets (the radio quiet sequence). In this scenario, there 
is still a universal correlation between jet power and radia- 
tive luminosity, but the same jet power corresponds to dif- 
ferent radio outputs in fast leptonic and slow mass-loaded 
jets. Possible observational tests of this hypothesis include 
measurements of circular polarization in the radio core, and 
of the location and intensity of the 511 keV annihilation line. 
For sources where the jet creates a plasma cocoon around the 
BH, comparison of the radio luminosity with the mechanical 
power necessary to inflate the bubble can also constrain the 
fraction of jet energy carried by protons and ions. 



far before stopping. A more likely explanation is that the 
lobes correspond to the locations where the heads of a pair 
of continuous jets (moving in opposite directions) encounter 
the interstellar medium and dissipate their bulk kinetic en- 
ergy. Our n ew radio observatio ns confirm the energy content 
inferred by lHardcastlel (|2005l ). that is ~ a few 10 erg in 
relativistic electrons. If the jet power is similar to the radia- 
tive core power, that is Pj ct ~ P ra d ~ 10 37 erg s" 1 (as is 
the case for example in Cyg X-1: Ru ssell et al J 120 07), such 
energy can be supplied in a few years. This would also imply 
that the lobe energy content, and therefore radio flux, can 
vary on a characteristic timescale of a few years, in response 
to changes in the injection power. 

However, this estimate is very uncertain because we do 
not yet know the fraction of kinetic and thermal energy 
stored in protons, ions and non-relativistic electrons. In sev- 
eral extra-Galactic jet sources, an independent estimate of 
the integrated jet power is available from the size of the 
jet-inflated cocoon around the BH. In those cases, the to- 
tal jet power can be up to a few hundred times higher than 
the po wer transferred to the radio-emitting relativistic elec- 
trons (|Cavagnolo et alj|20ld : ISoria et alj|2010l ; iPakull et ail 
I2010T I. Moreover, the radio core of GRS 1758-258 is 50 to 
100 times fainter, at a given X-ray luminosity, than the 
"standard" radio - loud class of Galactic BHs (Section 7.3; 
ISoleri et al] I2OI0I ; ISoleri fc Fenderll201ll ). We do not know 
whether that means that the total kinetic jet power (includ- 
ing both leptons and baryons) is correspondingly lower, or 
only the fraction of kinetic power carried into the lobes by 
the relativistic electrons is lower. But in either case, the ki- 
netic power energizing the radio-emitting electrons in the 
lobes is more likely to be < 10 35 erg s~ x . corresponding to a 
timescale > 10 3 yrs, compar able to the flux d ecay timescale 
due to adiabatic expansion l|Hardcastlell2005l ). 



7.4 Extended radio lobes 

Our 2008-2009 ATCA observations confirm the presence 
of extended radio lobes at apparently the same position 
as in 1992 and 1997, within the uncertainty due to dif- 
ferent beam sizes and sensitivities. Such structures are 
common in radio galaxies but much rarer a nd fainter in 
microquasars. It was suggested l|Heinzl [2051 ) that micro- 
quasars lack persistent lobes because they are mostly lo- 
cated in a low-density environment compared with AGN 
(scaled to their respective jet powers), so that any lobes 
formed have short lifetimes against adiabatic losses. Ejec- 
tions of fast-moving synchrotron-emitting blobs are well doc- 
umented in transient Galac tic BHs such as G RS 1915+105 
jRodrfguez fc Mirabel|[l999l ). GRO J1655-4 dTingav et al] 
19951). and XTE J15 5 0-564 l|Corbel et al1l2002t ): see "also 
Mirabel fc Rodriguez! (1999) for a comprehensive review. 
Such ejecta tend to decelerate and fade over shorter 
timescales (months or few years), after travell ing a char- 
acteristic distance « 0.5 pc (|Hao fc Zhand l2009). 



The distance between the core and the lobes of GRS 
1758—258 is an order of magnitude larger than that. Thus, 
they are unlikely to be individual blob ejections from a his- 
torical (unrecorded) outburst, unless the source is located in 
a very large cavity, that allowed the ejections to move that 



8 CONCLUSIONS 

We have studied the Galactic BH GRS 1758-258 in the 
X-ray and radio bands, during its spectral state transitions 
of 2000-2002, and during the hard state of 2008-2009. We 
showed that a compact jet is present during the hard state, 
and switches off in the soft state. We identified the X-ray 
states from their spectral and timing properties, constrained 
the disk parameters during the soft and intermediate state, 
and estimated a BH mass ~ IOMq and an X-ray luminosity 
« 0.01-0.03L E dd- 

We discussed its unusual X-ray state behaviour: the ap- 
parent luminosities in the hard and soft states overlap; in 
fact, it is often higher in the hard state. Since radiative effi- 
ciency in the hard state must be at least slightly lower than 
in the soft (thermal) state, this implies that the mass accre- 
tion rate in the hard state is often higher than in the soft 
state (assuming there is no beaming). This suggests that 
another physical parameter is more important than changes 
in the total mass accretion rate, for driving state transi- 
tions. The sequence of spectral variability in 2001 also sug- 
gests that a full disk was already present at least for several 
weeks before the main soft transition, when the X-ray source 
was in a hard or intermediate state. We conclude that the 
hard state of GRS 1758—258 may not be the same as the 
low/hard state of other BH transients successfully modelled 
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with truncated disks plus ADAFs; it may be an example of 
a Comptonizing corona/outflow on top of the disk. 

We showed that the radio lobes, first seen in 1992, were 
still present in 2009, apparently at the same location. This 
supports their interpretation as long-lived relic structures, 
perhaps associated with the interaction of the jet with the 
interstellar medium. We are leaving a more detailed study of 
the lobe evolution to a planned long-term ATCA radio mon- 
itoring of this source over the next few years. We measured 
the radio core luminosity in a series of unpublished VLA 
monitoring observations over 2001-2003, and compared it 
with the radio core fluxes measured in 1992-1993, 1997 and 
2008-2009. When the compact jet is detected, the 5-GHz 
flux typically varies between «s 0.1-0.5 mjy. Comparing this 
with the X-ray luminosity, we showed that GRS 1758—258 
belongs to the radio-faint sequence of Galactic BHs in the 
radio/X-ray plane. The dichotomy between radio-loud and 
radio-faint stellar-mass BHs for the same X-ray luminosity 
and spectral state is reminiscent of a similar dichotomy be- 
tween Seyfert and radio galaxies in the fundamental plane. 
BH spin is often suggested as a possible explanation, with 
fast-spinning BHs able to launch more powerful jets at the 
same accretion rate. We suggested an alternative explana- 
tion: both the radio-loud and radio quiet sequences may have 
the same jet power at a given accretion rate, but radio-loud 
BHs have a relativistic electron/positron jet, and radio-quiet 
BHs (such as GRS 1758—258) a sub- relativistic, mass-loaded 
jet, which is less efficient at emitting synchrotron radiation 
for the same kinetic power. 
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